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ABSTRACT 

^^ , Large scale correlations in the orientations of galaxies can result from alignments in their angular 

^^ ' momentum vectors. These alignments arise from the tidal torques exerted on neighboring proto-galaxies 

by the smoothly varying shear field. We compute the predicted amplitude of such ellipticity correlations 
using the Zel'dovich approximation for a realistic distribution of galaxy shapes. Weak gravitational 
Qh' lensing can also induce ellipticity correlations since the images of neighboring galaxies will be distorted 

^ , coherently. On comparing these two effects that induce shape correlations, we find that for current weak 

lensing surveys with a median redshift of Zm = 1, the intrinsic signal is of order 1-10% of the measured 
signal. However, for shallower surveys with z™ < 0.3, the intrinsic correlations dominate over the lensing 
signal. The distortions induced by lensing are curl-free, whereas those resulting from intrinsic alignments 
are not. This difference can be used to disentangle these two sources of ellipticity correlations. 

> 

^ I 1. INTRODUCTION 

Gravitational lensing can be used to map the detailed distribution of matter in the Universe over a range of scales 
Qv , (Gunn 1967). Systematic distortions in the shapes and orientations of high redshift background galaxies are induced by 
f~^ mass inhomogeneities along the line of sight. In strong lensing, a single massive foreground cluster will cause background 
f^ \ galaxies to be significantly magnified and distorted. Weak lensing, on the other hand, measures the cumulative effect of 
C^ • less massive systems along the line of sight statistically (Gunn 1967; Blandford et al. 1991; Miralda-Escude 1991; Kaiser 
1992; see a recent review by Bartelmann & Schneider 1999). 

The lensing effect depends only on the projected surface mass density and is independent of the luminosity or the 

dynamical state of the mass distribution. Thus, this technique can potentially provide invaluable constraints on the 

O ' distribution of matter in the Universe and the underlying cosmological model (Bernardeau, van Waerbeke & Mellier 1997; 

van Waerbeke, Bernardeau & Mellier 1999). There has been considerable progress in theoretical calculations of the effects 
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^ ' of weak lensing by large-scale structure, both analytically and using ray-tracing through cosmological N-body simulations 
'P. \ (Kaiser 1992; Bernardeau, van Waerbeke & MeUier 1997; Jain & Seljak 1997; Jain, Seljak & White 2000). 

Recently, several teams have reported observational detections of 'cosmic shear' - weak lensing on scales ranging from 

L^ ■ an arc-minute to ten arc-minutes (Van Waerbeke et al. 2000; Bacon, Refregier & Ellis 2000; Wittman et al. 2000; Kaiser, 

r\ ' Wilson & Luppino 2000). At present, these studies are limited by observational effects, such as shot noise due to the finite 

jrt number of galaxies and the accuracy with which shapes can actually be measured given the optics and seeing (Kaiser 

■ ■ ' 1995; Bartelmann & Schneider 1999; Kuijken 1999). In addition, the intrinsic ellipticity distribution of galaxies and their 

redshift distribution is still somewhat uncertain. These observational difficulties can be potentially overcome with more 

data. 

However, an important theoretical issue remains. In modeling the distortion produced by lensing, it is assumed that 
the a priori intrinsic correlations in the shapes and orientations of background galaxies are negligible. Correlations in the 
intrinsic ellipticities of neighboring galaxies are expected to arise from the galaxy formation process, for example as a 
consequence of correlations between the angular momenta of galaxies when they assemble. We compute the strength of 
these correlations in linear theory, in the context of Gaussian initial fluctuations. 

To do so, we approximate the projected shape of a galaxy on the sky by an ellipsoid with semi-axes a, b (a > b). The 
orientation of the ellipsoid depends on the angle ip between the major axis and the chosen coordinate system, while its 
magnitude is given by |e| — {a^ — b'^)/{a^ + fo^). Both the magnitude of the ellipticity and its orientation can be concisely 
described by the complex quantity e--°\ 



(°) = l^(°)|p2«V ^ r^(°) , ,-J°'l 



e^ ' = \C-'\e-'-^ = e 



+ 



(1) 



where the superscript '^°^ denotes the observed shape. 

In the linear regime and under the assumption of weak lensing, the lensing equation can be written as, 



e 



(o) = ' + 9 



where g is the complex shear and e the intrinsic shape of the source (Kochanek 1990; Miralda-Escude 1991). Furthermore, 
in the weak regime, correlations of this distortion field are 

(e(°)(xi)e(°)*(x2)) ^ (6(xi)e*(x2)) + {g* {^^)e{^^) + g{^^)e* {^^)) + (g(x03*(x2)) (3) 
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2 Spin Induced Galaxy Alignments 

where the * denotes complex conjugation.n In this paper we will examine the first term, which arises from intrinsic shape 
correlations. Previous analyses have focused on the third term of this expression, correlations due to weak lensing. The 
second term, which is due to correlations between the lensing galaxies and the intrinsic shapes of the galaxies being lensed, 
will not be addressed here. Naively however, we expect this contribution to be small, since the mean distance between 
the lensing and lensed galaxies far exceeds the distance scale over which angular momentum correlations are important. 
We will assume that shape correlations arise primarily from correlations in the direction of the angular momentum 
vectors of neighboring galaxies. Spiral galaxies are disk- like with the angular momentum vector perpendicular to the 
plane of the disk, so that angular momentum couplings will be translated into shape correlations. We will assume that for 
ellipticals the angular momentum vector also lies along its shortest axis on average, as it does for the spirals. However, 
since elliptical galaxies are intrinsically more round, the correlation amplitude will be smaller. Below we will use the 
observed ellipticity distributions of each morphological type in the computation of the shape correlations. For weak 
lensing, in contrast, the induced shape correlations are independent of the original shapes of the lensed galaxies. In 
the next sub-sections we will briefly review the origin of angular momentum and recent work on understanding intrinsic 
ellipticity correlations. 

1.1. Origins of Angular Momentum 
The angular momentum of the matter contained in a volume V is defined as. 



L(t) = / [r(i) - r(i)] X v(r, t) p(r, t) d'r, (4) 

Jv 

where f is the center of mass and p{r,t) is the density. Hoyle (1949) suggested that the origin of galactic angular 
momentum is tidal torquing between the proto-galaxy and the surrounding matter distribution. Most of the angular 
momentum of an object is imparted before the over-dense region completely collapses. After collapse, tidal torquing will 
be inefficient and the object will simply conserve its spin. 

Peebles (1969) used perturbation theory to calculate the growth rate of angular momentum contained within a comoving 
spherical region. For such a spherical region, there are no torques initially, so the growth occurs at second order as a 
result of convective effects on the bounding surface. In contrast, Doroshkevich (1970) showed that the angular momentum 
of a proto-galaxy grows at first order since, in general, proto-galactic regions are not spherical, generating an initial tidal 
torque. White (1984) described this process using the Zel'dovich (1970) approximation and showed that the spin grows 
linearly in time, for an Einstein-de Sitter universe. 

Following White (1984), we consider the growth of fluctuations in an expanding Friedmann Universe filled with pressure- 
free dust (p = 0) in Lagrangian perturbation theory. The trajectory of a dust particle can be written in comoving 
coordinates x = r/a in terms of the gradient of the gravitational potential '5, x(q, t) = q — Z?(t)V* (Zcldovich 1970). 
D(t) describes the growth of modes in linear theory and is proportional to the cosmological expansion factor, a(t), for an 
Einstein de-Sitter model. In terms of the Lagrangian coordinates q, the expression for angular momentum becomes 

L(t) = pofl^ / [x - x] X X d^g ~ poa^ / [q - q] x x (fq, (5) 

where po is the present mean matter density and Vl is the Lagrangian volume that corresponds to V . The latter expression 
is correct to second order since x — ~D{t)V'^ is parallel to the displacement. 

We can progress by expanding the gradient of the gravitational potential in a Taylor series around the center of mass, 

a„*(q) ~ a„*(q) + (q - Ci)pT^p (6) 

where the shear tensor is defined as the second derivative of the gravitational potential, Tai3{q) = da dp ^(q). The angular 
momentum of a collapsing proto-galactic region before turnaround then is given by, 

La = a^it) D{t) eafffTp^ la-y, (7) 

Icr-y is the moment of inertia of the matter in the collapsing volume, 

la^ = poa^{t) / (q - q)a(q - q)',d^q- (8) 

Note that the volume element that initially contains the matter is in fact much larger than that of the final galaxy in 
comoving coordinates. In this picture, the angular momentum is constant after turnaround. 

This formalism has been used to study how angular momentum arises during galaxy formation. Heavens and Peacock 
(1988) used Eulerian perturbation theory to compute the modulus of the angular momentum for galaxies, assuming the 
object to form at the peak of a Gaussian field (Bardeen et al. 1986). They found that there is a broad distribution in 
the angular momenta of collapsed objects, which is only weakly correlated with the heights of the density peaks around 
which galaxies form. 

Catelan & Theuns (1996a [CT96a]) expanded on this, working in Lagrangian space instead of Eulerian space. The 
results from these two approaches are very similar, but the resulting expressions are simpler for the Lagrangian case. 
CT96a approximated the shape of the object in Lagrangian space by an ellipsoid, which allowed the study of how angular 
momentum was correlated with other aspects of the matter distribution, such as its mass or its prolateness. The results of 

^ In the remainder of the text we will write two-point correlations functions in the following form; {f g') = (/(xi) (;(x2)) . 
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this analysis allows one to compute joint probability distributions, for example between the mass and spin of a halo. These 
were found to be in good agreement with the results from numerical simulations (Sugerman, Summers and Kamionkowski 
2000). 

Extending their approach, Catelan and Theuns (1996b) used second-order Lagrangian perturbation theory to estimate 
the contribution of non-linear effects, which they showed to be small. They also investigated the consequences of non- 
Gaussian primordial perturbations (Catelan & Theuns 1997), which they showed could have a significant effect on galactic 
spins. 

Lee & Pen (2000) re-examined the origin of angular momentum on galaxy scales and studied the statistics of both the 
magnitude and the direction of the present day spin distribution using numerical simulations. They developed a method to 
reconstruct the gravitational field, using only the direction of the angular momenta, since the predictions for its magnitude 
have a large variance. A central issue in determining the magnitude of L is the degree of correlation between the principal 
axes of the inertia tensor and gravitational shear tensor. CT96a attempted to take this into account around peaks, and 
found that such a correlation reduces the angular momentum by a small factor. Lee & Pen (2000) demonstrated using 
numerical simulations that this factor is in fact non-negligible. They however conclude that the approximation made by 
CT96a is adequate for determining the direction of the angular momentum vector but not for the magnitude. Here we 
extend the treatment of Lee & Pen (2000) and much of our notation and formalism follows their paper. 

1.2. Intrinsic Ellipticity Correlations 

There have been a number of preprints on this subject recently. We briefly review some of the results obtained by other 
groups here and we will compare our calculations and results with these in more detail in subsequent sections. 

Two groups. Heavens et al. (2000) and Croft & Mctzler (2000), have attempted to measure the strength of intrinsic 
correlations from high resolution cosmological N-body simulations (that evolve only the dark matter component) of the 
Virgo collaboration (Jenkins et al. 1998; Thomas et al. 1998; Pearce et al. 1999). Some assumption must be made 
to relate the dark matter halos in numerical simulations to the expected ellipticity of the luminous galaxies that form 
within them. Croft & Metzler (2000) measure the projected ellipticities of dark matter halos and the correlation of pairs 
as a function of separation. They then assumed that halo shapes are synonymous with galaxy shapes, and having done 
so claim to find a positive signal for the correlation on scales of the order of 20 h~^ Mpc (limited by the largest box 
size available). The results obtained in three dimensions were then projected into two dimensional angular ellipticity 
correlation functions, taking into account the viewing angle. They compute the induced correlations in the ellipticity and 
compare to recent reported measurements of the observed lensing signal. While there is a large uncertainty arising due to 
the unknown redshift distribution of the sheared background galaxies, they find that at most 10 - 20% of the measured 
signal could be attributed to contamination from residual intrinsic correlations. 

Heavens et al. (2000) have studied correlations in the intrinsic shapes of spiral galaxies also using the Virgo simulations. 
However, they use the angular momentum of the halo (rather than the actual shape as done by Croft & Metzler) and 
assume that its direction is perpendicular to that of a thin disk. They compute the 3-D ellipticity correlation function 
and its 2-D projection directly from simulations populated by ^ 10^ halos. They also conclude on comparing with recent 
measurements of the shear induced by lensing on large scales that the contamination from intrinsic correlations is small 
on most angular scales of interest - the contamination is roughly at the 10-20% level on scales of 0.1 — 10'. 

Catelan, Kamionkowski & Blandford (2000) have recently presented an analytic calculation to assess the importance 
of intrinsic galaxy shape alignments and the consequent mimicking of the signal produced by weak gravitational lensing. 
They make the Ansatz that the ellipticity is linearly proportional to the tidal shear and calculate correlations due to 
intrinsic shape correlations as a function of scale. (While originally meant to apply to ellipticity correlations resulting 
from angular momentum couplings, Catelan et al. now use this Ansatz only for ellipticities induced by the halo shapes. 
[M. Kamionkowski, private communication.]) They also consider possible means of discriminating the lensing signal from 
intrinsic allignments. 

Very recently, the first observational detection of the magnitude of spin-spin correlations has been reported by Pen, Lee 
& Seljak (2000). They construct the simplest quadratic two-point spin-spin correlation function in the context of linear 
perturbation theory and compare the statistic computed for galaxies in the TuUy catalog. They claim a detection at the 
97% confidence level out to a few Mpc. 

Several authors have pointed out that one of the important discriminants between the correlations arising due to lensing 
versus those from intrinsic alignments is the prediction of the existence of non-zero 'B-type' curl modes in the shear field 
in the intrinsic case (Kaiser 1992; Stebbins 1996; Kamionkowski et al. 1998). A detailed decomposition of the shear field 
into the 'B' and 'E', or pure gradient, modes for intrinsic correlations is presented in Crittenden et al. (2000). 

1.3. Schematic Outline 

Our goal is to calculate the two point correlation of the intrinsic shape distribution of galaxies, (ee'), as a function of 
projected distance. Since the following calculation is quite complex, we present a brief schematic outline to guide the 
reader and to clarify the simplifying assumptions that we make. Our approach is primarily analytic, but we also use 
numerical realizations of Gaussian fields to verify some of our results. 

The intrinsic ellipticity of a galaxy depends on its three dimensional shape, its orientation and on the direction of its 
angular momentum, e = e(S, L). Here we use S to denote the shape and orientation degrees of freedom. We will implicitly 
assume that the galaxy is ellipsoidal and that its angular momentum lies parallel to the shortest axis of the ellipsoid. The 



4 Spin Induced Galaxy Alignments 

expected correlation between ellipticities at different points is 

(ee'*)= /"dSdS'dLdL'e(S,L)e(S',L')-P(S,L,S',L'), (9) 

where V denotes the joint probability distribution. The present three-dimensional shapes of galaxies, quantified via 
their axis ratios, are primarily determined by 'local' processes like the extent of dissipation within the collapsing dark 
matter halo. Thus, we will assume that they are uncorrclatcd between neighboring galaxies, so that we can rewrite, 
7'(S,L, S',L') = ^(L, L')7^(S)'P(S'). For each galaxy, we can then integrate over all possible shapes and orientations 
to find the average ellipticity of a galaxy with angular momentum in a given direction, e(L) — J (iSe(S, L)7'(S). This 
integration is described in detail in Section 2. 
The resulting correlation is then simply given by, 

(ee'*) = I dt dt' e{t)e{t')V{U L')- (10) 

To proceed we need to understand the correlations between the directions of the angular momentum vectors of galaxies, 
or explicitly the nature of 7^(L,L'). Rather than attempt to calculate the angular momentum correlations directly, we 
instead relate them to correlations in the shear tensor, T, which yields itself more easily to linear theory. 

As discussed above, the angular momentum of a given galaxy depends on the tidal field and the moment of inertia, I, 
of all the matter that has turned around and which will eventually collapse to form the galaxy. To compute 7^(1/, L'), one 
needs the full joint probability function, 7^(1, T, I', T'). We make the simplifying assumption that the moment of inertia at 
a given point is significantly correlated only with the shear at that point, so that 'P(I, T, I', T') = 'P(I|T)'P(I'|T')'P(T, T'). 
Given a form for 7'(I|T), one can derive 7'(L|T). However, since the local stress tensor depends on the details of the 
mass distribution outside the collapsed object as well, 7^(I|T) is not accurately known. In Section 3, we follow LPOO and 
assume that 7^(L|T) is Gaussian, and use the most general form that the correlation matrix could have as a function of 
the shear tensor. This allows us to derive an expression for e(T) = J (iLe(L)7'(L|T), the expected mean ellipticity for a 
given shear tensor. 

With these assumptions, the ellipticity correlation depends only on how the tidal field is correlated from place to place. 



= /"dTdTe(T)e(T')7'(T,T') =F(C(r,,)), (11) 



where C{vij) is the correlation matrix of the shear tensor, which will be Gaussian distributed if the underlying fluctuations 
are Gaussian. The ellipticity correlation is now only a function of the separation r^j. In Section 4, we calculate the 
correlations of the shear tensor as well as the moments required to find (ee'). Later in Section 4, we also examine how the 
correlations of the shear change if they are sampled only at peaks of the density. This is to account for the fact that we 
are sampling galaxies, which do not form at random positions in space. 

Till now, the ellipticity correlations we have been considering are in three dimensions. These correlations must be 
projected into two dimensions to compare with weak lensing predictions and measurements. In Section 5, we do this 
projection using Limber's equation (Limber 1953). This allows us to take into account the clustering of galaxies whose 
ellipticities are sampled. 

In Section 6, we examine the implications of our results for weak lensing observations and the prospects for measuring 
the intrinsic signal in on-going surveys like the SDSS and 2dF. We conclude in Section 7, with a more detailed discussion 
of our assumptions and the uncertainties involved in our calculations. 

2. INTRINSIC SHAPE DISTRIBUTIONS 

In this section, we relate the observed, projected shapes of galaxies to their three-dimensional shapes in the absence 
of lensing. When calculating intrinsic shape correlations, it is important to take into account the distribution of three 
dimensional shapes since the strength of the signal depends strongly on it. For example, a spherical galaxy will appear 
round when viewed from any angle, consequently its presence will tend to suppress intrinsic shape correlations. We first 
consider the simplest case, where the galaxies are modeled as thin disks, with the angular momentum vector perpendicular 
to the disk plane. This is a fairly good approximation for spiral galaxies. We then consider the effects of projecting more 
realistic galaxy shapes, modeling them to be tri-axial with Gaussian distributed axis ratios. 

For a galaxy with a thin disk, the exact dependence of the observed ellipticity is easy to calculate. The shape of 
such a disk-like galaxy depends strongly on the observing angle, appearing round (e — 0) when viewed face on and very 
elongated (e — 1) when viewed edge on. When viewed from an angle with respect to the perpendicular, the disk will be 
foreshortened by a factor of cos in one direction. The magnitude of the observed ellipticity is then, 

, , 1 - cos2 1 - L2 

I'l ^ TT^ ^ TTIf' "'> 

where L = (sin 9 cos (^, sin 9 sin 0, cos 9) , is the unit spin vector. The observed ellipse has its long axis oriented perpendicular 
to the projected angular momentum vector, so that tj) — (p + 7r/2. 

For realistic galaxies, however, the relation between the observed and intrinsic shapes can be much more complicated. 
The finite thickness of the galaxy puts an upper limit on how elongated the observed shape can be. Observed galaxy 
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Fig. 1. — Left panel: The average ellipticities of the different morphological types of galaxies seen from a given angle 9 with respect to 
the angular momentum vector. Also shown is the result for a sample with all types weighted by the observed fractions in the APM survey. 
All scale roughly as the thin disk case, Right panel: The distribution of ellipticities for the different morphological types based on the LBL 
intrinsic shape distributions. The differences between the LBL model and the fit from lensing surveys (dashed curve) could either be due to 
there being a different morphological mix in high redshift surveys or evolution in the intrinsic shape distributions, particularly for the spirals. 



samples also contain a mix of morphological types - ellipticals, spirals and spheroidals - each of which has a different 
distribution of intrinsic shapes. 

We consider the intrinsic shape distributions found by Lambas, Maddox & Loveday (1992; LBL hereafter), extracted 
from the Bright Galaxy Survey of the APM catalog. They used triaxial models to describe the observed ellipticity 
distributions of the various morphological classes and obtained fits for the distribution of the underlying axes ratios. They 
assumed Gaussian distributions for the scaled axes ratios {a ~ \) for all three morphological classes of the form. 



V{b) (X exp 



2al 



V{c) ex exp 



2a2 



(13) 



truncated such that < c < 6 < 1. For the spiral population, they found it necessary to include the effects of a finite disk 
thickness in order to explain the deficit in the high ellipticity tail of the observed distribution. The best-fit parameters 
for spirals were found to be 6o = 1-0, Cb = 0.13, cq = 0.25 and Cc — 0.12. LML also demonstrated that simple oblate 
or prolate models were not capable of reproducing the observations for elliptical galaxies and that triaxial models were 
required. For ellipticals, the best-fit parameters were found to be &o = 0.95, cri, = 0.35, cq = 0.55 and CTc = 0.2. Finally, 
the best-fit parameters for spheroidals were found to be Jdq = 1.0, ai, = 0.3, cq = 0.59 and CTc = 0.24. 

Stark (1977) derived the relation between the three-dimensional axes ratios and the ellipticity, which we adapt to the 
case at hand. Knowing the distribution of galaxy shapes for a given galaxy type, we can calculate the average ellipticity 
of a galaxy with angular momentum at an angle 6 with respect to the line of sight. 



e(L)- / dS,e(S„L,)7'(S,) = 



d0 / dbV{h) I dcV{c)e{9,(t),b,c), 
Jo Jo 



(14) 



where e{9, (j), b, c) is complex and given in Appendix A. We have performed these integrations numerically and display 
the results for the amplitude in Figure 1. By symmetry, the average orientation angle is perpendicular to the projected 
angular momentum. 

As can be seen in the left panel of Figure 1, e is different for each morphological type. The maxima, corresponding to 
when the galaxies are seen edge on, are significantly less than 1 due to the finite thicknesses of the galaxies. Interestingly, 
however, the dependence of e on 6* for each of the morphological classes roughly scales identically to the disk case. That 
is, it is a good approximation to assume the same functional dependence on angle with an overall scaling: 



HL.)\ 



a- 



1+Lr 



(15) 



where < a < 1 is a measure of the relative galaxy thickness. For the more realistic distributions, a ranges from 0.85 for 
spirals, down to around 0.5 for spheroids and ellipticals. 

The mean redshift of the APM survey is Zm = 0.1 and the composition of the sample is roughly 10% ellipticals, 25% 
spheroidals and 65% spirals. While these fractions might be representative of a local field sample, the redshift distribution 
of background galaxies of interest in lensing studies is considerably higher, so the morphological mix could be much 
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different. To examine if this is the case, we plot both the distribution of ellipticities for the LBL populations and that 
measured for Icnsing studies 

T'del) oc |e|e-(l^l/"-3'"" (16) 

(Brainerd et al. 1996; Ebbels et al. 2000), which is represented by the dashed curve on the right panel. The LBL 
distributions provide a poor fit to the results of the higher redshift field surveys. This could be due to two effects: either 
a different morphological mix or a different, perhaps non-Gaussian, distribution of intrinsic axis ratios for spirals. In 
addition, since the higher redshift surveys are more likely to be dominated by irregulars, the intrinsic shape distributions 
are expected to evolve. 

Though we find that the LBL distributions on the whole provide a poor-fit to that from Icnsing studies, the inferred 
mean value for |e| is in good agreement. For a thin disk, lejrp = 7r/2— 1, so our simplified model has |e| = Q;|e|TD = 0.57a. 



The mean of the measured ellipticity distribution (Eqn. 16) is 0.42, implying that a = 0.73. This is consistent with the 
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mean value for a. computed from the APM sample on the left panel of Figure 1. In the following sections, we will be 
explicitly computing correlations between the components of L in order to calculate (ee') based on the definition in Eqn. 

(15). 

3. ELLIPTICITY AND THE TIDAL FIELD 

We wish to relate the ellipticity directly to the tidal field. The easiest way to do so is to consider the real and imaginary 
pieces of the distortion field separately: 

I cos(2<^) ^ « : ;:;^ ,;: = c£^^ 

1 + Ll 

I sin(20) = a ^~r^ .^"^! - a-^ . (17) 

1 + LlLl + Ll 1 + Ll 

When the observation frame coincides with the frame where the stress tensor is diagonal, (LyLx) — and the expected 
distortion is purely real. Next we need to relate the angular momentum to the tidal field. 

One central assumption in this work is that the expectation value of the angular momentum at a point is solely a 
function of the tidal tensor at that point. As the next step in calculating ellipticity correlations, we need to understand 
this relationship more quantitatively. That is, we need to know the probability of a given spin direction for a specified 
shear field, or effectively the form of 7'(L|T). As discussed above, the angular momentum of a collapsing region is given 
by La — a?{t) D{t) ea/s-y Tf}„ I^-y ■ The crucial issue is how the moment of inertia for a collapsing object is related to the 
tidal field that it experiences. Unfortunately, this requires a precise understanding of what determines the region which 
eventually collapses into the galaxy, which in turn depends on the positions of nearby over-densities. This remains a 
major unsolved problem. 

Catelan and Theuns (CT96a), studying the variance of the amplitude of the angular momentum of galaxies, initially 
assumed that the tidal tensor and the moment of inertia are entirely uncorrelated. However, since galaxies form preferen- 
tially at density peaks, CT96a also consider the suppression of angular momentum that arises around peaks in a Gaussian 
due to correlations between the inertia and the shear. In both of these cases, if one considers the frame where the inertia 
tensor is diagonal, the off-diagonal terms of the shear tensor are expected to be uncorrelated with the inertia tensor. 
However, around peaks the amplitude of these off-diagonal terms is suppressed, which results in lower angular momenta. 

As we will show, the amplitude of the angular momentum has little effect on the magnitude of ellipticity correlations, 
which is chiefly determined by correlations in the directions of the spins. For simplicity, we will assume each component 
of the inertia tensor to be Gaussian distributed, so that the resulting distribution for the angular momentum given some 
shear tensor is also Gaussian distributed, 

P(L|T) = ^^e-^-e^-^"/^ (18) 

(27r)2|Q|2 

where Qaf3 = (LaLp) is the correlation matrix. In the frame where the shear is diagonal, the inertia tensor is uncorrelated 
with it and the correlation matrix has the form, 

"(T33-T22)2 

{LaLp)<x (Tii-r33)2 . (19) 

{T22~Tn)\ 

This has the same form in the peaks case studied by CT96a. The relation can be rewritten as 

(LaLp) = (L^) (ha(3 - fa^f^A , (20) 

where T is the unit normalized traceless tidal tensor {TapTaf} = 1). Note that the angular momentum is independent 
of the trace of the tidal tensor, so we can consider T to be traceless without loss of generality. Our final expression for 
the ellipticity correlations will be independent of the proportionality factor (L^), so uncertainties in this factor will be 
irrelevant here. 
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It has been argued that the approximations made by CT96a underestimate the correlations between the moment of 
inertia and the tidal field, and result in overestimating the angular momentum that is produced when compared to 
simulations. Lee and Pen (LPOO) suggest that the CT96 approximations may do well in predicting the direction of the 
angular momentum, but not its amplitude. They consider the most general correlation between the shear and inertia 
tensors, 

(LaL/s) = {L^) f — :r—5ap - aTa^f^p j , (21) 

where < a < 1. The CT96a case of uncorrelated moment of inertia corresponds to a = 1. In the extreme a = case, the 
direction of the angular momentum is random, independent of the tidal tensor. Ironically, stronger correlations between 
the moment of inertia and the tidal field make the direction of the expected angular momentum more, not less, random. 
The directions can be further perturbed by non-linear interactions, particularly if the magnitude of the angular momentum 
is small originally. LPOO investigate this in N-body simulations and find that the relation is best fit by a = 0.24. 

The distribution of the direction of the angular momentum vector is given by integrating over the amplitude of the 
momentum, 

V{t\'T) = [ L^dLViL\T) (22) 



CO 



L^dL- ^ 



(2^)*|Q|5 

^ ALcQ-fM-i, (23) 



4^121 

where Q = Q/{L'^). The variance of the two point expectation value of the direction of the angular momentum is then 

1 

■47rTQ| 
As shown by LPOO and in Appendix B, for small a this implies that 



L^Lp) ^ I d^L L^Lp^^-^{La,Q-'pLp)-2. (24) 



(LaL/s) == o (1 r)*^"/? + -rTa-,T.yf). (25) 

Combining the above results, we can compute the average ellipticity for a given shear tensor at one point, 

1{T)= j dte{t)V{t\T) 

= aJd^L ^^^^J^^ [cos20 + zsin2^]^-l^(£„Q-^L;3)-i (26) 

For small a, we can approximate Q^^ ~ ^[Sa/s ~ a(<5a/3 ~ STa-yTyp)]. Therefore, {LaQ'^lLf))^^ ~ |Q|2[1 + i^ — 
^LaL/^Ta^T^/j]. Inserting this in the integral, by symmetry the surviving terms are 

e(T) = — -^ / d9 / d(t>{cos 20 + i sin 20) Ifi^f^i cos^ + f2jf^2 sin^ + 2fi^f^2 sin cos ( 

Stt Jq 1 + cos-^ Jq I 



9/*7r '5/1 /*2'7r 

aa / ,„ sm 



/ d9 — / d<?!)(cos2(/) + isin20) [A + Bcos2(/) + Csin20] (27) 

Stt Jq 1 + cos-^ Jq 

where A = \{fi^f^i + T27A2), B = \{fi-,f-,i - f2-y%2) and C = fi^f-,2. Finally, using ^^'^ dO jfj^ = tt - 8/3, we 
find that the integral of Eqn. (26) evaluates to 

Qtt 1 . . . . . . 

e(T)-e+(T)+zex(T)-aa(6-— )[-(ri^T^i-r2^T^2) +^Tl^T^2]. (28) 

The numerical factor, ^ — 6, is very nearly unity and we shall drop it for convenience here. Thus, the average ellipticity 
for a given tidal field is quadratic in the tidal field and is suppressed both by a factor due to the finite thicknesses of the 
galaxies (a) and by the randomization of the angular momentum vector (a) . 

4. CORRELATIONS IN THE TIDAL FIELD 

In the previous section we related the ellipticity to the tidal field, and here we calculate the moments of the tidal field 
required to derive ellipticity correlations. Since the tidal tensor is the second derivative of the gravitational potential, its 
statistical properties are directly related to those of the matter density. We assume that the underlying density field is 
Gaussian distributed, which implies that the tidal field is also Gaussian. However, the unit normalized tidal field which 
is of relevance to us will not be Gaussian distributed. 

Knowing the full statistics of the tidal field, it is possible to calculate expectation values of observables such as the 
ellipticity correlation (ee'*) = / (iT(iT'e(T)e*(T')7'(T, T'). Since the ellipticity is quadratic in T, we need to compute 
linear and quadratic two point functions of the normalized shear field, (TT') and (TTT'T') respectively. We first compute 
the correlations of T which are necessary to evaluate these moments. 
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4.1. Correlations o/T 

To begin, we compute the two point expectation value of the tidal tensor in terms of the power spectrum of fluctuations. 
A Fourier expansion of the shear tensor yields, 

Ta^(x) = -J d^k K ki3 ^(k) e^'^"', (29) 

where \l/(k) is the Fourier transform of the gravitational potential, which has a power spectrum defined as, 

(vl/(k) *(k')) = fe(k - k') P^(fc). (30) 

^From this, it is straight forward to calculate the two-point correlation function of the tidal tensor, 

C„;5^,(r) = (r,^(x)T^,(x')> = /c^3kfc„fc^fc^fc,P^(fc)e*''•^ (31) 



where the separation is r = x — x' Jj 

To evaluate the correlation function, it is useful to relate Fourier space components back to the real space derivatives 
via ika = da- Thus, we have 

= 2TTdadpd^da I dkk^P^{k)jo{kr). (32) 



Here we have performed the integration over the angular directions of the Fourier modes, and jo{kr) = sin kr/kr is the 
zeroth order spherical Bessel function. 

It is useful to rewrite the derivatives as da — {dr/dxa){d/dr) — XaD, where the operator D = {l/r){d/dr). Using the 
identity daXp = Sap, we find 

Cai3ja{r) = 27ra„a/3 f dkk^P^{k){S^^D[jo{kr)] + x^x^D^Mkr)]) 

= 2Tr[SapSja + SajS/Ba + SaaS/Jj] / dk k'^ P^ {k) D'^[jo{kr)] 

r 

+ 27r[rarfj6jcr + Tar^dpa + rar^Sp^ + rpr~^5aa + rpr^Sa^ + r^r^dafi] I dk k^ P^ik) D'^[jaikr)] 

+ 2Trrarpr^ra / dkk^ P^{k) D'^[ja{kr)]. (33) 



It is possible to use Poisson's equation to substitute the power spectrum of the potential with that of the density, k^P^i(k) = 
Psik), in units where AirGpo — 1 and G is the gravitational constant. Using the identity D" jo{r) = (—1)" r~^ jnir), the 
above simplifies to 

Ca/37cr(r) = [6a/3S.ycr + Sa-y^lSa + 5aa5f3-f]C2{r) + raffir^r^ (iir) 

+ [fa-rpS-^a + raf^Spa- + faf„5p^ + fpf^Saa + fpf^Sa^ + f^f^Saf}] CsC*^), (34) 

where 

C„(r) = (-1)"^ fdkk--^Ukr)Psik), (35) 



r 



and jn{kr) is the n*'* spherical Bessel function. This is identical to the expression derived in LPOO, wherein it was shown 
that these C, functions are related to the density correlation function, ^(r) — J dk k^ Ps{k)jo{kr) , and integrals of it. 

If the density field is smoothed, its correlation function levels off as r ^ 0. In this limit, (3 and (^4^-0 and the above 
expression reduces to 

Ca0ja{O) = C2{0)[Saf3Sycr + SajS/3a + SaaS/3j]- (36) 

Since C2(0) — ^(0)/15, this corresponds precisely to the variance of the tidal field found by CT96a (equation (38) in 
Appendix A.) The correlation function simplifies dramatically when averaging over directions f, 

1 /■ 1 

This is useful when the correlation length is much smaller than the depth of the survey. 
^Note that the positional vectors are in Lagrangian space and were denoted by q in Section 1. 
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4.2. Correlations o/ T 

In this sub-section, we calculate the two and four point moments of the unit normalized traceless tidal field, T. Since 
T is not Gaussian distributed, these moments are not necessarily simply related. In the next sub-section, we use these 
results to derive the final ellipticity correlation. 

For simplicity of notation, it is useful to treat the six degrees of freedom in the shear tensor as a vector, which we 
shall denote with capital Roman subscripts, T = (Tn, T22, 233, T12, T13, T23). While using this notation, it is important to 
remember that the shear transforms as a tensor under rotations, rather than as a vector. To further simplify the notation, 
we shall write the shear field at a displacement r from the origin as T'. Thus in this notation, the correlation matrix 
becomes a six by six matrix, (TaTq) = [Cr]AB- 

Since the tidal field is Gaussian, with a two-point correlation matrix given by C, we can write the expectation value 
for an observable like TaTb as. 



{TaT's) = 



rTATUe—- 



1 „-if'^c~^f 



(38) 



(39) 



|C|l/2(2^)6- 

where T = (T, T'). The matrix C has the following block diagonal form, 

„ _ Co Cr 

^ — /"( /"( 

Or Oo 

where Co is the zero-lag correlation matrix and Cr contains the two point correlations. For galaxies that are separated 
by distances greater than the smoothing scale, we can assume that Cr ^ Co and expand in powers of C^Cq to invert 
C to second order: 

Cg (1 -|- CrCg CrCg ) ^'"'O '"'r'-'o 

— Cq CrCg Cq (1 -|- CrCg GrOg ) 

We perform a Taylor expansion of the exponential. 



c-i 



(40) 



exp[-if ^C-if] ~ [1 + TCo"^ Cr Co' T' + ^{T C^' C, C^' T'f 



1 



cxp[--(TCg-^r + r'Cg-^r')]. (4i) 



To evaluate the linear two point function of T, we must keep terms to first order in CrCg . The expectation value can 
then be written as, 

{TaTb) = 



Ta f'g [T C^ ^ CrCg" ^ T'] e- ^ (^ ^o ' ^+^' ^o ' ^') 



Co|(2^)6^"^«^^ ^" ^''^° 

{fATc)[C^'CrC^']cD{f'sTD') 



where 



{Ta Tb) 



(fT 



f^TBe-*(^^o-^^), 



(42) 



(43) 



(27r)3|Co|i/2- 

is proportional to the mean magnitude of the tidal field. This is evaluated in Appendix C by transforming variables, 
Ta ^ Raa'Ta', to a basis in which the correlation function is proportional to the unit matrix. Using the results derived 
there, the linear two point function is shown to be, 

<^-^-'>-225^[^r]- (^^) 

where Cr is the correlation function of the traceless part of the tidal field. Though this was evaluated in the large 
separation limit, its value at zero lag is very close to the exact result. 

The quadratic two point function of T is evaluated in an analogous way, except here we must keep terms to second 
order in C^C^ . Making this substitution, the quadratic two point function is 



{TaTbT'cT'o) 



(PTdPT' 

|C|i/2(27r)6 



-' rf^l „-iT" C^'T 



TATBTcToe—- 



= {Ta Tb) {T^T'b) + ^{Ta TbTeTf) {fc TdTgTh) [C^^ C, C^']eg[C^' C, C^']fh. (45) 



As before, the expectation value of 

{TaTbTcTd) 



(fT 



Y/jTa TbTcTd e 



-l(TC-^T) 



(46) 



(27r)3|Co|i/2 

is computed in the transformed basis and is derived in Appendix C. The final form of the quadratic two point function is 

then 

2 



{TaTbT'cT'b) = 



1 



[Cr]Ac[Cr]Br> + [Cr] AD [Cr] BC 



local terms. 



(47) 



,14C2(0), 

where Cr is defined as above. The local terms correspond to the reducible parts of this fourth order moment and do not 
contribute to the ellipticity correlation. 
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Fig. 2. — Left panel: The computed three dimensional ellipticity correlation function averaged over angles and plotted as a function of 
separation. The signal is appreciable at separations smaller than the smoothing scale (< l_Mpc) and falls off as ^^(r). For comparison we 
have plotted the exact derived function [solid curve] against the approximation of Equation fell) [dashed curve] . Right panel: The expectation 
value of the two components of the two point ellipticity correlation (e-|-e+) and (ex^x) divided by the approximate curve (solid curve of left 
panel) for various values of the viewing angle 9 = 0, 7r/4 and it/2 are shown. Note that the 2 components are equal for 9 = 0. 



4.3. Correlations of the Ellipticity 

Using the results derived earlier, we are finally in a position to calculate ellipticity correlations. Recall that the ellipticity 
correlations are given by 



ee 



dtdt'e(t)e* (t')7'(t, t') 



-{Ti-fT-fi 



T2^f^2) + ifi^f^2][^{f[„T'^i -T2aTa2) 



^T[X2])- 



(48) 



(e+^V) 



We choose the separation vector r to lie in the x — z plane at an angle 9 from the line-of-sight which we assume to be 
parallel to the z— axis. This choice implies that frames in which the ellipticities are measured lie parallel to the projected 
separation. 

Inserting equation ( [47| ) the two non-zero components of the ellipticity correlation are, 

W(i4^)' [336C|(r) +472C2(r)C3M + 155Cf (r) + 58C2M C4(r) + 26C3(r) C4(r) + 3C|W 

+ 4(l8C2(r)C3(r)-7C|(r)-8C3WC4M- CIM) cos 2 

+ {nC.l{r)+&Ur)Ur)+&Ur)U{r)+Cl{r)) cos40], (49) 

(exe'x) = 4#(T4c^)M42C2'(r) + 59C2(r) C3(r) + ISClW + 5C2(r) C4(r) + C3(r) C4(r) 

+ (9C2(r)C3(r)+5C|(r) + 3C2(r)C4(r)- C3(r)C4(r))cos20]. (50) 

For a simple model with £^{r) (x l/r, (e+e+) and (exe'x) £^rc plotted in Figure p^, computed assuming top-hat smoothing 
on a 1 ft,~^ ]Mpc scale. 

These functions are explicitly anisotropic and depend on the angle between r and the line of sight. Much of the angular 
dependence can be understood intuitively by considering the symmetries of the problem. When the line of sight is parallel 
to r (6* = 0), there is no longer a distinction between e_|_ and ex, so that {t+f'+) and (cxCx) are identical. In addition 

-9 and ^ tt — so that the only surviving terms are 



and 



these correlations are invariant under the transformations 9 
either constant or proportional to cos26' or cos 46*. This anisotropy is demonstrated in the right panel of Figure ]Li 
Asymptotically at large r, the angle averaged behavior is approximately 

^" '~ 84 ^2(0)- 



(51) 



This is very close to the exact expression at large r, as is shown in the left panel of Figure 4.3, The factor of 1/84 has been 
derived for large separations. At zero lag, it can be computed directly from the fourth moment of the unit normalized 
traceless shear tensor and is found to be 1/60. 
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Fig. 3. — The angle averaged ellipticity correlations from Gaussian field realizations on a 512"^ grid compared to the analytic results. The 
correlations at peaks of the density field (shown as solid squares) appear to match well with the results from random positions (the dot-dashed 
line). At zero-lag they both asymptote to the exact result (1/60) shown by the horizontal dot-dashed line. The full analytic correlations 
(solid curves) as well as the analytic estimate where the large scale power has been removed (dashed curve) to account for the finite size of 
the realizations are both shown. These are not valid below the smoothing scale. The analytic results match well in the region where they are 
valid, i. e. at scales larger than the smoothing length, but note the importance of large-scale power (here we have set a = a = 1). 



4.4. Peaks in a Gaussian field 

Galaxies do not form at random positions, but at peaks of the density field (Bardeen et al. 1986). It is possible that 
this sampling could bias the expected correlation of galaxy ellipticities since our analytic correlations have been computed 
for random points. We examine such a potential bias using numerical realizations of Gaussian fields, and also use these 
to check the validity of our analytic results. 

We create realizations of Gaussian fields on a 512'^ grid with a power spectrum Pg oc fc^^, corresponding to a density 
correlation which falls off as 1/r. The density field is smoothed with a spherical top-hat filter of approximately four grid 
units. Peaks are identified as positions where the density field exceeds the value at each of its six nearest neighbors. The 
tidal field is calculated at each point using differencing. At each point, we compute and subtract the trace and finally 
unit normalize the tidal tensor. 

We checked first that the moments of the normalized tidal field match our analytic expectations at zero lag. For 
example, we can analytically calculate the variance of one component of the tidal tensor, (T^j,). From isotropy, this can 
be shown to be 2/15, which we have verified in the realizations. We have also checked numerically other exact quadratic 

and quartic relations such as (f^^) = 1/10, (f^^) = 4/105 and {f^^f^y) = 1/105. 

Using the definition in Equation (48), we compute the correlation function of the (angle averaged) ellipticity on the grid 
both for peaks and random field points, assuming a = a = 1. These are compared with the derived analytic results in 
Figure 3. The correlation function for the peaks (^ 70,000 in the box) and the field points are in excellent agreement for 
both small and large separations. At zero lag, they both asymptote to the analytic result of 1/60 (marked in Figure 3 by 
the horizontal long-dashed line). The maximum deviation occurs around the smoothing scale and is of the order of 10%. 
This implies that while peaks might preferentially be sites of galaxy formation, as far as intrinsic ellipticity correlations 
are concerned, there is no substantial bias between peaks and random field points. 

The numerical correlation function starts dropping below the exact analytic expectation (plotted as the solid curve) for 
large separations. This is due to missing large-scale power on the grid: an analytic calculation that incorporates the same 
lack of power on large scales as the grid is shown as the dashed line. Its agreement with the numerical results demonstrates 
that the steeper fall off of the numerical estimates is indeed an artifact of the finite box size. Clearly this is a worry for all 
numerical computations of the ellipticity correlation function. Additionally, on separations smaller than the smoothing 
scale, our analytic estimates fall below the numerical results and the exact values derived from the one point moments. 
This discrepancy arises because the approximation we made in computing the analytic results, C(r)/C(0) ^ 1, is invalid 
on small scales. 

5. THE PROJECTED CORRELATION FUNCTIONS 

Up to now we have been focusing on the three dimensional ellipticity correlation function. While this is in principle 
observable (Pen et al. 2000), weak lensing studies usually consider the ellipticity correlation projected onto the sky. The 
projection of the intrinsic signal will enable us to compare directly with weak lensing estimates and judge its importance 
as a possible contaminant for these measurements. 
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5.1. Limber's equation 

We begin by making some general comments about the basis dependence of the two dimensional correlation functions. 
Previously, we calculated the three dimensional correlation functions in a special basis, one where the separation vector 
was coplanar with one of the axes used to define the ellipticity. We showed that the ellipticity correlations in this basis 
were functions only of the 3-D separation and the angle cos 9 ^ z ■ r. In two dimensions, this basis is equivalent to taking 
one axis vector to be parallel to the 2-D separation f. n In this special basis, which hereafter will be denoted by the 
superscript r, the ellipticity correlations will only be functions of the distance between the two points. We will denote 
these functions as, 

e+(|r1) ^ {el{x)elix + r)); CAM) ^ {e^-A^VA^ + r^) ■ (52) 

The cross correlation, (e!j_(a;)ex(^+ ^))j is zero due to parity (x — > —x) invariance. 

The basis for these correlation functions depends on the separation vector, and thus on which pair of galaxies one is 
considering. It is often useful to work in a fixed basis on the sky for the ellipticities. The ellipticity measured in an 
arbitrary basis with an angle (j) relative to the separation vector is given by: 

e+ = e'j_ cos 2(j) — e^ sin 20; 

ex = €+ sin 2(/) + e^ cos 2(j>. (53) 

In such a fixed basis, the correlation function depends on cf) as, 

Ci(|fU) = (e+eV> =e+(H)cos2 20 + ^AM) sin^ 2cf> 

C2(|r1,</>) = (exe'x) =ex(|r1)cos^20 + f+(|r1) sin^ 20, (54) 

for all pairs separated by |r| at an angle of (p with respect to the chosen basis. The sum of these is a function of the 
separation only, while the difference has a simple dependence on 0: 

Ci(|rU)+C2(|r1,0)=C+(k1)+ex(|r1); Ci(|rU) - Csdrl, 0) = K+(|r1) - ^xdrD] cos40. (55) 

Recent measurements of the shear from weak lensing have focused on the variance of the magnitude of the ellipticity 
averaged over a patch, which depends only on the sum. 

We next consider the projection into two dimensions and use an approach similar to that used by Heavens et al. (2000) 
and Croft & Metzler (2000). Assuming that we are working on a small area of the sky, the observed patch of sky is 
approximated by a plane. The projection uses Limber's equation to take into account the clustering of galaxies, 

/ zfz^dzidz2ijizi) ip{z2) [1 + £,gg{r)] (e+(xi)e+(x2)) 



MM) - 



/ zfz^dzidZ2ll^izi) ^P{Z2) [1 + Cgg{r)] 



p n^s ^ I zlzldzidz2ip{zi) 'IP{Z2) [1 + Cggjr)] (cx (xQex (xa)) 

^'^' '^ Jzfzidz^dz2i^{z,)i^{z2)[l + Cgg{r)] ' '^^^> 

where |r|^ = (r^ — (zi — 22)^), ip{z) is the observational selection function and (,gg{r) is the galaxy-galaxy correlation 
function. 

Note that while the ellipticity correlation function is calculated in Lagrangian coordinates, the projection is performed 
in Eulcrian space. The Eulerian separations will differ from the Lagrangian ones due to peculiar velocities of galaxies, 
which we have ignored here. For galaxies near each other, the Eulerian separations will in general be smaller than the 
corresponding Lagrangian ones. This will result in a suppression of the intrinsic projected ellipticity correlations. We 
expect this effect to be small at large separations, but it could be significant closer in. 

5.2. Qualitative Features 

In Figure wc plot the projected correlation function calculated for the simple model in which the density correlation 
falls off as l/r. We have assumed a top-hat smoothing scale of 1 h~^Mpc. The clustering term is taken to be of the 
form ^gg = [r/Rs)~^ where Rs = 5h~^ Mpc is the clustering scale and P = 1.8 (Loveday et al. 1992). Finally, following 
Heavens et al., the selection function is taken to be ip{z) = ^-i'^l'^o) which has a mean redshift of z„i = 1.4zo. In the 
figure we plot the functions ^+ and ^x for a = a = \. They are nearly identical at small angular scales, but deviate at 
larger separation. 

The qualitative features of these correlation functions can be understood by looking at the various scales in the problem: 
i?m, the mean depth of the survey; i?^, the clustering length of galaxies and i?Oj the smoothing scale. If the density 
correlation function falls off as l/r, the three dimensional ellipticity correlation scales as ^_|_ ~ ^x — a^a^i?o/84(i?Q 



where we have used the approximation from Eqn. (pl|). Both the numerator and the denominator of Equation (56) 
contain a clustering term, which dominates at small angular separations and a mean contribution. 

The denominator in Limber's equation is essentially the average number of galaxies within an angular distance Q from 
a given galaxy out to the volume of the survey. At large separations, the denominator scales as the volume squared, or 

i?^. Corrections from clustering are of order R^{Rs/Rm)^d^~^ , which become important at angles 9 < [Rs/Rm)^^ ■ 
For a survey with median redshift 2;,„ ~ 1, this occurs at about 2 arc seconds. For a shallower survey with Zm ^ 0.1, this 
occurs at much larger scales, of order a few arc minutes. 

^For clarity, we use the bold face for three dimensional vectors and the notation r for 2-D vectors. 
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Fig. 4. — The projected 2-d ellipticity correlation functions, §+ a^nd ^x i ior various values of the median redshift of the distribution of 
galaxies. Here we have assumed the galaxies are well described by thin disks and are perfectly aligned with the linear predictions (a = a = 1.) 
but deviate at larger separation. Both increase strongly with decreasing redshift, roughly proportional to z^ . 



The numerator in Limber's equation is the projected cUipticity correlation function, weighted by the number of pairs 
at a given separation. Again, the clustering term dominates at small scales. For separations less than 9 '^ Ro/Rm, the 
three dimensional correlations are effectively constant with an amplitude of 1/60, so the behavior is identical to that of 
the denominator, a'^a^R^{Rs/Rm)^9^~^ /60. At very large separations, 6 > Rs/R,n, clustering is not important but the 
numerator also falls off inversely with angular separation, RqR^O^^ /84. Between these regimes, Ro/Rm < & < Rs/Rm, 
there is a transition where the clustering contribution falls off quickly. 

Thus the projected correlation functions have a number of distinct regimes. At very small separations, clustering 

dominates both the numerator and the denominator, leaving the correlation constant (a^a^/60 for 9 < {Rs/ Rm)^~^ ■) 
On slightly larger scales, but smaller than Ro/Rm, clustering dominates the numerator, but not the denominator, and 
the correlation falls off as a power law, ^^ ^ a^a^ /6O{Rs/Rm)^0^~^- There is then a brief transition region where the 
correlation falls fairly quickly. Finally on very large scales, the mean values dominate both the numerator and the 
denominator and the correlation falls off as a^a^0~^i?§/(84i?^). 

It is straight forward to understand the dependence of the correlation functions on the mean redshift of the survey. The 
typical 3-d separation of galaxies with a given angular separation is directly proportional to the survey depth i?,„. Thus 
if the three dimensional ellipticity correlations fall off as r"", then the projected correlations fall off as i?~". For the case 
we have been considering, ^ ex 1/r, so that the ellipticity correlations fall off as l/r^, and the projected correlation drops 
as z~^. This is clearly seen in Figure ^. 

6. INTRINSIC ALIGNMENTS VERSUS WEAK LENSING 

In the previous sections we presented an analytic expression for the intrinsic ellipticity correlation function, which we 
now evaluate for realistic surveys. We also compare this with recent measurements of cosmic shear and with theoretical 
weak lensing predictions at high and low redshifts. 

The amplitude of intrinsic correlations depends on both the mean thickness of galaxies and on their degree of alignment 
with the tidal field. In Section 2, we argued that observed shapes of galaxies are characterized by a = 0.73. The degree 
of alignment of the galaxies with the predictions from linear theory, parameterized by a, was measured by Lee and Pen 
(2000) in N-body simulations and found to be fairly small, a ~ 0.24. This implies that the correlations could be suppressed 
by non-linear effects. 

However, we are interested in the correlations of spins with each other, not necessarily in how they align with the 
predictions from linear theory. The LPOO measurement was a one-point measurement, and thus can not account for 
'correlated randomizations.' Non-linear interactions between galaxies could lessen the correspondence of their spins with 
the linear predictions without changing how well the spins correlate with each other. Thus using the one point value 
of a will likely underestimate the amplitude of the spin correlations. The measurements of the three dimensional spin 
correlations from the Virgo simulations (Heavens et al. 2000) can be used to measure an effective a which takes this effect 
into account. At 1 Mpc, they find a correlation of approximately 5 x 10^^ which is in remarkable agreement with that 
found by Pen et al. observationally. When compared to the analytic prediction of a^a^/60, this yields an effective value 
of a = 0.55. (Heavens et al. treat the galaxies as thin disks, so that a = 1.) Here we will present results for both a = 0.24 
and a = 0.55 to demonstrate the possible uncertainty of our predictions. 

In the left panel of Figure H, we plot the sum of the intrinsic correlation functions, C+ + Cx , for a median redshift of 
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Fig. 5. — The intrinsic correlation signal versus the predictions from weak lensing and current observations. Left panel: i;+{6) +Cx{^) 
for a median redshift of 1, compared to the measured shear correlation function. At small separations, the intrinsic signal is approximately 
one percent of the measured value. The amplitude depends on the value of the assumed average galaxy thickness (a) and the parameter a 
that describes how well the angular momentum of the galaxy is correlated with the shear field. We plot a = 0.24 (full line) and a = 0.55 
(short-dashed line) which correspond to the values inferred from numerical simulations by LPOO and Heavens et al. (2000) respectively. 
a = 0.73 corresponds to the value determined from the observed distribution of ellipticities (Ebbels et al. 2000). The data are: van Waerbeke 
et al. (2000) - solid squares ; Wittman et al. (2000) - filled circles ; Kaiser et al. (2000) - open circles; and Bacon et al. (2000) — filled 
triangle. The long-dashed line is the theoretical prediction from Jain & Seljak (1997) computed for a Q\ = 0.7 galaxy cluster normalized flat 
universe, ~ 4.75 X 10~^(6'/arcmin)~''*^. Right panel: as in the left panel but for the predictions for a shallower survey such as SDSS and 
2dF with median redshift Zm = 0.1. The intrinsic signal is again shown for two values of a, and the theoretical prediction for weak lensing is 
the long-dashed line (for Zm = 0.1) and dotted-long-dashed (for Zm = 0.5). The lensing prediction for Zm = 0.1 is extrapolated from the Jain 
& Seljak fit beyond the stated range of validity. For such low rcdshifts, the intrinsic signal is significant and may dominate over the lensing 
contribution for most scales. 



1, a galaxy smoothing scale of lh~^ Mpc and the parameter choices described above. We also show the measured shear 
variance compiled from the recent literature. The weak lensing prediction from Jain & Seljak (1997) for an f^A = 0.7 
galaxy cluster normalized flat universe fits the data well. At small separations, the intrinsic signal contaminates the 
lensing one at the level of a few per cent, modulo the uncertainties in a and a. 

Note that there can be ambiguities in plotting measurements of the cosmic shear. One issue is whether the correlation 
function or the tophat variance is plotted. For a simple 1/6 correlation function, the variance is nearly a factor of two 
larger than the correlation at the same scale. In addition, some authors plot the variance as a function of the tophat 
smoothing radius, while others instead plot it as a function of the diameter. Finally, some authors quote the variance 
of each component of the complex shear field, while others quote the variance of the modulus of e. Here we plot the 
theoretical predictions for ^_|_ + ^x, the correlation of the modulus of the ellipticity. In contrast, when plotting the data, 
we have used the modulus variance for a given tophat radius. Since the correlation function implied by the data is slightly 
lower than the variance, the relative contribution of the intrinsic correlations is somewhat larger than is naively implied 
by the figure. 

The relative importance of the intrinsic correlations increases dramatically as the depth of the survey is reduced. As 
the observed galaxies are closer to us, the lensing signal falls because there is less intervening matter to lens them, while 
the intrinsic signal grows since the galaxies are physically closer to each other for a given angular separation. Jain & 
Seljak (1997) show that the lensing signal scales as z^^"^ for a flat SIa = 0.7 universe. In contrast, if the density correlation 
function is proportional to r", then the ellipticity correlation scales as z^". On galaxy scales, n ~ — 1, hence the intrinsic 
amplitude grows rapidly since the signal scales as ~ z~^. 

In the right panel of Figure |5| we compare the intrinsic correlations for a shallow survey with Zm = 0.1 such as SDSS or 
2dF to the lensing signal expected from the theoretical analysis of Jain & Seljak (1997). We show their fits for Zm = 0.5 
and also extrapolate the fit to Zm = 0.1. The latter is beyond the stated range of validity, but should give an approximate 
idea of the lensing amplitude. For low Zm < 0.3, the intrinsic signal is significant and may dominate over the lensing 
contribution on most scales. Clearly, large surveys like 2dF and SDSS offer exciting possibilities for measuring intrinsic 
shape correlations. 

There are other important distinctions between the lensing and the intrinsic correlation signals. For example, the lensing 
signal depends on the amplitude of the mass fluctuations, parameterized by as- In contrast, the intrinsic correlations 
depend only on correlations of the direction of the shear field and are therefore largely independent of the amplitude of 
the fluctuations. 

In addition, another difference arises in how the intrinsic signal depends on morphological type. Weak lensing is in some 
sense democratic, as all galaxy types are distorted in the same way. This is not the case for intrinsic correlations however. 
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We have shown that this signal depends on the distribution of axis ratios. Spiral galaxies are characterized by a ~ 0.85, 
while ellipticals and spheroidal galaxies typically have a ~ 0.5. In addition, the alignment of the angular momentum with 
the shortest axis is likely to have more scatter in elliptical galaxies, resulting in an effective lowering of the value of a. 
Thus we expect the intrinsic correlation to be suppressed by more than a factor of two. This hypothesis can be checked 
observationally by using color criteria to separate the morphological types of galaxies, since ellipticals tend to be redder 
than spirals. 

7. SUMMARY 

In this paper, we have presented a calculation of intrinsic correlations in the observed ellipticities of galaxies resulting 
from angular momentum couplings. We have focused on the angular momentum which arises in linear theory and is 
associated with the local tidal field. The three dimensional spin correlations were projected using Limber's equation to 
obtain the expected 2-d ellipticity correlations. These intrinsic correlations were shown to dominate over the weak lensing 
signal for shallower surveys. 

A number of assumptions were made in order to make the calculation tractable. Foremost of these is the assumption 
that angular momentum plays the central role in aligning the observed galaxy shapes. Other factors, such as the initial 
distribution of matter which fell in to form the galaxy, could also conceivably have contributed to the observed shapes. 
However, the angular momentum is special in that it is approximately constant during the later evolution of the galaxies. 
Galaxies typically have had many dynamical times to virialize, and we expect most of the dependence on the initial matter 
distribution to be lost. This is particularly true for spirals, and holds for ellipticals which are slow rotators and have spin 
parameters of the order of 10%. Although their rotational time scale {^ 1 Gyr) is much longer than their dynamical 
time (~ 100 Myr), they have undergone enough rotations in a Hubble time to erase any memory of the alignment of the 
principle axis (Dubinski 1992). 

At small separations, other factors, such as the recent history of galaxy collisions, might also affect the ellipticity 
correlations. In addition, it is important to remember that we are probing only the light distributions, which reflect the 
matter distribution only at the very central parts of the galaxies. 

It is essential to understand precisely how the ellipticity correlations depend on the angular momenta. The dominant 
contribution to the correlations comes from alignments in the orientations of the galaxy ellipticities. The elongations of 
the galaxy light distributions are expected to be orthogonal to the direction of their projected angular momenta. The 
magnitude of the ellipticity may to some extent depend on the magnitude of L; for example, galaxies with larger angular 
momenta may appear more disk-like. Even so, the form of this relation is largely irrelevant for understanding the ellipticity 
correlations. This is because the galaxy orientations are expected to be isotropic on average, so there must be correlations 
in the alignments for ellipticity correlations to occur. However, the magnitude of the ellipticity has a significant mean 
value. Therefore, galaxy alignments already introduce ellipticity correlations even in the absence of correlations in the 
magnitude of the ellipticities. 

To see this, consider the ellipticity correlation 

(6e'*) ^ (|,||,'|)(e^^('^-'^')) = [|ef + {{\e\ \-e\){W\ ~ mKe''^'^-'^'^)- (57) 

The first relation follows from assuming that the magnitudes of the ellipticities are independent of their orientation 
correlations. Recall that the distribution for ellipticities described in Eqn. (jlfl) has a large mean value, |e| ~ 0.42. The 
variance of this distribution is significantly smaller than the square of the mean, ((|e| — |e|)^) = 0.055 ~ 0.3|ep, so that 
even perfect correlations between magnitude of the ellipticities would only result in a small modulation of the overall 
correlation. 

Since the ellipticity is proportional to e^"'^, it is quadratic in the angular momentum components perpendicular to the 
line of sight (Eqn. (|l7|).) The ellipticity correlation is therefore quartic in the angular momenta: (ee') oc {LLL'L'). This 
result should be contrasted with the Ansatz of Catelan, Kamionkowski and Blandford, which assumes the correlation to 
be quadratic in the angular momentaJj 

The correlation strength (ee') also depends on the mean ellipticity, which in turn depends on the galaxy type. Spiral 
galaxies are more flattened than elliptical galaxies, and thus will have a larger correlation. We have also assumed that 
the angular momentum is parallel to the shortest axis of the galaxies, which should be a good approximation for spirals, 
but may not be as good for elliptical galaxies and could suppress their correlation further. 

Another major simplifying assumption we have made is that linear theory is sufficient to calculate these angular 
momentum correlations. We might hope that this is a good approximation, since most of the angular momentum is 
expected to be imparted before the object starts to collapse and enters the non-linear regime. While there are non- 
linear corrections, N-body simulations have shown the linear approximation to be surprisingly robust (Lee and Pen 2000, 
Sugerman et al. 2000). We have attempted to account for the effects of non-linear evolution by parameterizing the extent 
to which spin alignments are suppressed in comparison with linear predictions. This parameter was estimated from the 
N-body simulations of Lee and Pen (2000) and Heavens et al. (2000), as well as by comparison to measurements of 
observed ellipticity correlations seen in the TuUy catalog (Pen et al. 2000). 

The amplitude and shape of the ellipticity correlation function can be understood intuitively. Recall that the ellipticity 
is a function of the shear tensor, which is the second derivative of the potential. By virtue of Poisson's equation, the trace 
of the shear tensor is the density. Therefore we expect the correlation of the other components of the shear field will drop 

*These authors have recently re-examined this issue and now find results consistent with those we have presented here [M. Kamionkowski, 
private communication.] 
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at the same rate as the density correlation function. Since the ellipticities are quadratic in the shear field, correlations 
in them will fall as the density correlation function squared, (ee') oc ^^. The order-of-magnitude of the amplitude of the 
correlation at zero-lag follows from simple symmetry arguments. The shear tensor has six degrees of freedom, but only 
five are relevant since angular momentum is independent of the trace. The typical magnitude of a fourth order moment of 
a unit tensor in a five dimensional space is 1/35. Therefore, these considerations suggest that the ellipticity variance will 
have a comparable amplitude. Fuller consideration shows that, from Equation(|5l]), (ee') oc 1/84 x ^'^ (r) / S,'^ {r = 0). Our 
analytic calculations are valid for random points in a Gaussian field, but galaxies are usually assumed to form at density 
peaks. We performed large realizations of Gaussian fields and checked that our results are good approximations for such 
special sampling. 

We have compared the strength of the intrinsic correlation to that expected for weak lensing. The intrinsic signal grows 
as the depth of the survey decreases, because then galaxies close on the sky are on average also physically closer together, 
hence they are more correlated. The weak lensing signal, on the other hand, becomes weaker, since there is less matter 
between us and the lensed objects. For surveys typical of weak lensing, with a median redshift of z^ — 1, the intrinsic 
signal is of order of 1 per cent of the weak lensing amplitude. However, for shallower surveys such as SDSS or 2dF, the 
intrinsic signal may dominate the lensing one, on small scales. Therefore, SDSS and 2dF are ideally suited for studying 
intrinsic correlations in the orientations of galaxies. 

The intrinsic ellipticity depends on the square of the tidal field, whereas the lensing distortion is linear in the shear. 
As a direct consequence, the distortion field is curl-free when induced by lensing, but not when intrinsic correlations are 
present as well (Crittenden et al. 2000). The detection of such 'magnetic' modes will be an invaluable way of separating 
lensing from intrinsic correlations. 

Finally, in this paper we have concentrated on intrinsic correlations of galaxies. Applying a similar reasoning to clusters, 
one could hope to study the shear field on much larger scales. The alignment of clusters of galaxies is dominated by the 
intrinsic alignment of the major shear axis. Their dynamical time is longer, and they form later, so we would expect the 
initial formation alignment to persist, implying ellipticities linearly proportional to the shear. The correlation should then 
drop as the correlation function ^(r) instead of its square as is the case for spin alignments. The qualitative features are 
reported for Gaussian random fields in Pen (2000) and for simulations by Tseng and Pen (2000). 

We thank L. van Waerbeke for useful conversations. RC and TT acknowledge PPARC for the award of an Advanced and 
a post-doctoral fellowship, respectively. PN acknowledges support from a Trinity College Research Fellowship. Research 
conducted in cooperation with Silicon Graphics/Cray Research utilizing the Origin 2000 supercomputer at the Department 
for Applied Mathematics and Theoretical Physics (DAMTP), Cambridge. 
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APPENDIX 



FROM INTRINSIC TO PROJECTED SHAPES 

Here we derive an expression for the projected ellipticity, e = jeje^*''', for a general ellipsoid when viewed from an 
arbitrary angle. We follow the treatment of Stark (1977) and consider a galaxy as an absorption-free stellar system, in 
which the volume brightness is constant on similar ellipsoids. 

The general equation for the ellipsoid of the constant volume brightness in the coordinate frame of the galaxy (x, y, z) 
is, 

t^ x^ + u^ y'^ + z"^ ^ al, (Al) 
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where t is the axis ratio c/a, u the axis ratio c/b and a^ is the variable that parameterizes the vohimc brightness. (Note 
that one axis ratio u used here differs from the one used by LBL, which is h/ a.) We want to transform the above equation 
into a frame that is ahgned with z along the line of sight, which is accomplished by a general rotation, characterized by 
the first two Euler angles and 0, 

- sin (h cos 9 sin 6 sin 9 
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Stark shows that projecting the volume brightness along the line of sight yields curves of constant surface brightness 
described by. 



al = {l)x'' + 2q)x'y' + {L)y'^ 

where as parameterizes the surface brightness and, 

/ EE /(0, 9, t, u) = t^ sin^ 9 sin^ </» + u^ sin^ 9 cos^ </) + cos^ 9 

j = j{(j),9,t,u) = t^u'^ sin^ 9 + f cos^ (j) cos^ 9 + u^ sin^ (j) cos^ 9 

k = k{(j), 9,t, u) = {u^ — t^) sin cos (/) cos 6' 

l = l{(j), 9, t, u) = t'^ sin^ + u^ cos^ 0. 

For a given set of axes ratios and observation angle, j, k, I and / are constant, so that the projection of curves with 
constant surface brightness are similar ellipses. Therefore the projected image of a galaxy which has luminosity constant 
on similar ellipsoids has isophotes which are similar ellipses, with the same position angle. 
These isophotes correspond to ellipses with 



(A2) 

(A3) 
(A4) 
(A5) 
(A6) 
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l + \/T^ 



^1 + e' 



where 7 = At^v?'f/{j + /)^, and /3 the ratio of the short axis to the long axis, and 



^ 



1 



2fc 



(A7) 



(A8) 



the angle between the major axis and the x' direction. Therefore the general expression for the projected ellipticity of a 



galaxy with axes ratios t and u seen from a line of sight 



with respect to the galaxy frame is, e(0, 0, i, u) = ^l — je'^'^'^ . 



MOMENTS OF L 
As discussed in the text, the variance of the expectation value of the direction of the angular momentum is 

{LLf3)= j cfLL^Lp-^{UQ-lLfi)-i. 
J 47r|Q|2 



(Bl) 



Here, we estimate it in the limit of a ^ 1. A similar discussion can be found in LPOO. Writing the measure <PL = dLid(j) 
the angular integral in the above expression can be evaluated explicitly. In the frame where Q is diagonal, 

dl, I d0_--^--^_^--^, (B2) 



(iiLi 



where 



B = Q22 cos^ (, 



4TT\Q\i iALl + B)i 



Q33 sin^ 



Substituting tan^ tp = ^L^, we have. 



(iiii) = 2 



The tp integral can be evaluated exactly to be 
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(B3) 



(B4) 



(B5) 



The approximation is valid in the small angle limit tan ^ y^A/B <C 1, which corresponds to assuming that a in cqn. (9) 
is small. Using tan^"*^ tp ~ ip — tp^ /S, the integral becomes, 

{L,Li)o,l r d<p-±-^[B-i ^AB-i], (B6) 
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since, 



we have, 
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To linear powers in a, the first term in the integral becomes, 
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Similarly, the second term in equation (B6) gives, 
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Thus, (LiLi) ~ i(l 

elements remain zero. Thus the full correlation matrix becomes 



^■^Tii). Similar expressions hold for the other diagonal correlations and the off-diagonal 
full correlation matrix becomes, 

. . 1 ^n 'An - 

(Bll) 



(1 - -rl'^ap 



3a ~ ~ 

O 



3' 5 

THE LINEAR AND QUADRATIC SHEAR TWO-POINT FUNCTIONS 

Here we perform integrations useful in evaluating the two and four point functions of T. As described in the text, we 
will transform to a new basis, T = ((Tn + T22 + T^i)lVl, (Tii - T22)/V2, (Tn + T22 - 2r33)/V6, V2T12, V^Tu, V2T23). 
This is a convenient basis to integrate over the trace (since it is irrelevant to the angular momentum) and because the 
correlation function has a particularly simple form. The transformation matrix between these bases is: 



^ ^ " 
— — — — 

V2 
V2 
V2 

We will work in this basis throughout this appendix. 

Let us consider first the correlation function at zero separation in this new basis. In terms of its original indices, 

Co = ^^(0)[(5q,^(5^ct + '^a-ySfia + ^acrSfij]- In thc ncw basis this becomes 



1 1 
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[Ca]AB — Raa'[Co]a'B'[R ]b'b 



^diag(5,2,2,2,2,2) 



(C2) 



The factor in the exponential of the Gaussian distribution, TC^^^T, can be written in this basis as TCq^T = ^^^(0)(TrT^ + 
15|Tp/2), where |Tp = J2^a=2 '^a i^ the modulus of the traceless part of T. 
In this basis, it is simple to calculate (T4 Tb) as 

{-TaTb)^ I ,^jZ ,,„ fArBe-h\r\o)iTrT^+',\r\V2)], (C3) 

TrTfJi/j/vo we can easily perform 



(27r)3|Co|i/2 
Converting thc measure to (fT = -^(rYT:T\T\'^d\T\d'^f and rewriting To = \T\fD 



the integrations. The trace integral yields 



while thc modulus integral is 



r|5d|r|e-5[i5C-^(o)|r|V2)] ^ 8(2e(0)/15)3. 



(C4) 
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The determinant in this basis is simply |Co|^^^ — 4\/l0(C(0)/15)'^ so that we find 

,^^, 2 /e(o)v /• ,4^^ ^ 2 (my^^\ 16 (my, .p.^ 

(^Ar^)=^(^^j ydrr^r^=-^(^— j —^^^ = ^^-^ (^_ j 5^^, (C6) 

where (5^15 runs only over the non-trace indices (2-6). Effectively, operating (T^Tb) on a vector projects out the trace 
part of the vector. We are finally in a position to find the linear two point function: 

64 64 

SAc[Cr]cDSDB = -, r ^rn\ \-^r]AB, (C7) 



~ 15^(0) "^^'^^^"""^ 15<(0) 
where the tilde denotes that the trace has been projected out of the correlation function. (If P is the projection operator, 
then C = PCP^ . In the original basis, this projection operator is R^^SabR) At zero separation, this gives an answer to 
within 10% of the exact value (T^Tb ) = 5ab/^, a- remarkable fact when one remembers that this was derived assuming 

Moving on, we next try to evaluate the quadratic two-point correlation in this basis. 



{MBfh%) = J ^^i^^AfB%%e-h^' 



= {fAfB){%fi,) + \{fAfBrErp) {fcfDrcrH) [C(^c,cy]EG[c^^c,cy]pH. (C8) 

Now, at zero-lag the quartic moment in the transformed basis, 

{fATBrcrn) = / (2,)t|J„|i/2 ^^^Brcr,,e-^K-^(")(^^^+^^l^l^/^)l. (C9) 

Again making the substitution Tc — TTc H — \=TiT 5ic, the surviving terms are of the following form, 

. J=dTrT\T\'^d\T\d^f TrT^ 

{TaTbTcTd) = j ^ (2^)3|c„|i/2 rATEW^TcTo + ^Sicho]- (CIO) 

As above we can perform the integrals simply. For the first term, the trace integral is identical to equation (C5), while 
the modulus integral is 

|r|6d|r|e-i[i5«"'(o)l^lV2] ^ i5^(2^(o)/i5)i. (Cll) 

Jo 2 

The angular integral yields 

/o 2 
d'^ffATBTcfo = —ttASab^cd + Sac'Sbd + SadSbc]- (C12) 

Again the indices range over 2-6, since the trace has been effectively projected out. For the second term, the angular 
integral is as it was for the linear two point function (C6), while the trace integral becomes 

dTrT(TrT)2e-i[«"W(TrT^)] ^ ^2^^(0)1, (C13) 

^ — OO 

and the modulus integral gives, 



|r|4d|r|e-^[i5«"'(")l^l'/2l ^ ^(2e(0)/15)i (C14) 

2 



Putting these all together, we find that 

^2 



{fAfBTcTD) =^ 

15 



„{Sab6cd + SacSbd + SadSbc) + SABSicho 



(C15) 



Finally, we are in a position to evaluate the quadratic two point function. This expression has a number of terms, but it 
can effectively be broken into a local part, which includes terms proportional to Sab or Scd, and a non-local part. Since 
only the latter terms contribute to the ellipticity correlation, we will keep only these here. This non-local part is a simple 
function of the correlation of the trace- free components, 

{fAfBfifD)=(Y^^j [CacCsd+CadCbc] + local terms. (C16) 

In the limit of small separations, the total (local and non-local) correlation function should approach [SabScd + Sac^bd + 
SadSbc]/S5. 



